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The Riemann - Silberstein - Majorana - Oppenheimer approach to the Maxwell 
electrodynamics in vacuum is investigated within the matrix formalism. The matrix 
form of electrodynamics includes three real 4x4 matrices: {—iBq + aWj)^{x) = 0, 
where ^{x) = (0,E(x) + icB(x)). Within the squaring procedure we construct four 
formal solutions of the Maxwell equations on the base of scalar Klein - Fock - Gordon 
solutions: {^f^, = {ido + a^di + G?di + a^^s) ^>(x), where ^>(x) satisfies 

equation d°'da^{x) = 0. The problem of separating physical electromagnetic waves 
in the linear space {Aq^'^ + \\^^ + \2^'^ + As^*^} is investigated, several particular 
cases, plane waves and cylindrical waves, are considered in detail. 

1 Introduction 

Special relativity arose from study of the Maxwell equations symmetry with respect to motion 
of references frames: Lorentz [Ij, Poincare j2j, Einstein ^ Naturally, an analysis of the Maxwell 
equations with respect to Lorentz transformations was the first objects of relativity theory: 
Minkowski [4J, Silberstein [5l[6], Marcolongo [8], Bateman [9], and Lanczos |10j . Gordon 
Mandel'stam - Tamm [T2l [T3l [Ti] . 

After Dirac [15j discovery of the relativistic equation for a particle with spin 1/2 much work 
was done to study spinor and vectors within the Lorentz group theory: Moglich [16], Ivanenko - 
Landau [E] , Neumann [18] , van der Waerden [19] , Juvet [20] . As was shown any quantity which 
transforms linearly under Lorentz transformations is a spinor. For that reason spinor quantities 
are considered as fundamental in quantum field theory and basic equations for such quantities 
should be written in a spinor form. A spinor formulation of Maxwell equations was studied by 
Laporte and Uhlenbeck [21], also see Rumer [28]. In 1931, Majorana [23] and Oppenheimer 
[22] proposed to consider the Maxwell theory of electromagnetism as the wave mechanics of the 
photon. They introduced a complex 3-vector wave function satisfying the massless Dirac-like 
equations. Before Majorana and Oppenheimer, the most crucial steps were made by Silberstein 
[5] , he showed the possibility to have formulated Maxwell equation in term of complex 3-vector 
entities. Silberstein in his second paper [6] writes that the complex form of Maxwell equations 
has been known before; he refers there to the second volume of the lecture notes on the differential 
equations of mathematical physics by B. Riemann that were edited and published by H. Weber 
in 1901 [7]. This not widely used fact is noted by Bialynicki-Birula [103j ). 

Maxwell equations in the matrix Dirac-like form considered during long time by many au- 
thors: Luis de Broglie [MUSS [301 |36j, Petiau Proca [271 E], Duffin [29j, Kemmer [3ll El E] , 
Bhabha [32], Behnfante [33l[34], Taub [35], Sakata - Taketani [37], Schrodinger [39l|40], Heitler 
[12], [la He], Mercier [17], Imaeda [H], Fujiwara [50], Ohmura [5l], Borgardt [SH [59], Fe- 
dorov [53] , Kuohsien [54] , Bludman [55] , Good [56] , Moses [571 EOl [76] , Lomont [58] , Bogush 
- Fedorov [Hj, Sachs - Schwebel [66], Ellis [68], Oliver [70], Beckers - Pirotte [71], Casanova 
[72], CarmeH [73], Bogush [74], Lord [75], Weingarten [77], Mignani - Recami - Baldo [78] . 
Newman [79], [80], [82], Edmonds [83], Silveira [86]; the interest to the Majorana - Oppen- 
heimer formulation of electrodynamics has grown in recent years: Jena - Naik - Pradhan |87j . 
Venuri [88j, Chow [89j, Fushchich - Nikitin [90j, Cook [92l|93], Giannetto Yepez, Brito - 



Vargas [27], Kidd - Ardini - Anton Recami [22], Krivsky - Simulik [M], Inagaki [TU2] . 
Bialynicki-Birula [l03l lOl 128] , Sipe [TOS], [l06], Esposito [I08], Dvoeglazov [IM] (see a big 
list of relevant references therein)- Gersten |lU7j . Kanatchikov [llOj . Gsponer |112j . Ivezic 
[TT31 [TM [TT51 [TTBl [TT71 [TTSl [TT91 [Ml [m], Donev - Tashkova. [I25l [T26l[T27] . 

Our treatment will be with a quite definite accent: the main attention is given to possibilities 
given by the matrix approach for explicit constructing electromagnetic solutions of the Maxwell 
equations. In vacuum case, the matrix form includes three real 4x4 matrices a^: 

{-ido + a^dj)'^{x) = 0, 

where ^{x) = (0, E(3;) + icB(x)). With the use of squaring procedure one may construct four 
formal solutions of the Maxwell equations on the base of scalar solution of the Klein - Fock - 
Gordon equation: 

{^0^ ^1, ^r^i ^ ^-Q^ ^ ^iQ^ ^ ^2g^ ^ ^3^^^ ^ = . 

The problem of separating physical electromagnetic solutions in the linear space {Ao^'^ + Ai^*^ + 
A2^2 + As^^} is investigated. Several particular cases are considered in detail. 



Complex matrix form of Maxwell theory in vacuum 

1129] (with the use of usual notation 



Let us start with Maxwell equations in vacuum [38 
for current 4- vector = (p, J/c) , = 1/eo/io )'■ 



div cB = , rot E = 



dcB 



div E 



rot cB 



dct ' 
(9E 

eo dct ' 



Let us introduce 3-dimensional complex vector iIj^ = + icB^ , with the help of which 
the above equations can be combined into (see Silberschtein 0(6], Bateman [9], Majorana [23] . 
Oppenheimer [2^, and many others) 



(2) 



let xq = ct, Oq = cdt- These four relations can be rewritten in a matrix form using a 4- 
dimensional column ip with one additional zero-element [90[ I130j : 



{-ido + a^dj)^ = J 







a 



10 
-10 

1 



a 



a 






1 











1 


-1 











-1 





{a' 


f = -I, 











a 



J= — 







1 

-1 



2\2 



a 



J 

if 
if 
1 
-1 


2^2 _ J 

2 



[a 



a 



(3) 
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3 Method to construct electromagnetic solutions from scalar 
ones 



The above matrix form of Maxwell theory : 

{-ido + a^dj)-^ = , 







(4) 



permits us to develop a simple method of finding solutions of Maxwell equations on the base of 
known solutions of the scalar massless equation by Klein - Fock - Gordon. Indeed, in virtue of 
the above commutative relations we have an operator identity 



{—ido + a^di + 0^82 + a^^s) {—ido — a^di — 0^82 



i-dl + df + dl + dl) 



Therefore, taking any special scalar solution (— c^o +^1 +^2 +^3) ^{x) = , one can immediately 
construct four solutions of the Maxwell equation: 

{-ido + a^di + a^d2 + 0^83) = , 

where are columns of the matrix 

{ido + a^di + a'^d2 + a^ds) <^{x) = 





^1$ 


82^ 


83'^ 






-83^ 


82<^ 


-92$ 




i8o^ 


-8i<^ 


-53$ 


-82^ 


81^ 


i8o^ 



Thus, we have four formal solutions of the free Maxwell equations (let Fa{x) = 8a^{x)): 



iFo 


Fi 


F2 


F3 


-Fi 


iFo 


-F3 


F2 


-F2 


F3 


iFo 


-Fi 


-F3 


-F2 


Fi 


iFo 



(6) 



4 Electromagnetic plane waves from scalar ones 

Let us specify this method for an elementary example, starting from a scalar plane wave prop- 
agating along the axis z: 



$ = Asm{ujt — kz) = Asm.{koXo — ^3X3) = Acos^p . 



(7) 



The recipe ([5]) gives 



{^"} = A 



iko —k3 

iko ^3 

—k3 iko 

/c3 iko 



cos If . 



(8) 
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Because the left and the right columns have non- vanishing zero-component, they cannot repre- 
sent any real solutions of the Maxwell equations. However, two remaining seem to be suitable 
ones (the factor A cos 93 is omitted) : 











iko 




+ icB^ 


-h 




+ icB^ 







E^ + ciB^ 



II 



A 











^3 




E^ + icB^ 


iko 




E'^ + icB^ 







E^ + ciB^ 



Thus, we have two wave solutions: 

/ E'^ = —k^A cos 99 , 

II E^ = k^A cos (f , 
In both cases, the ratio E/B equals to the speed of light: 

h 



B^ 



B' 



ko 
c 

ko 
c 



A cos ip , 
A cos 93 



(9) 



E 
B 



ko/c 



Besides, both waves have are clockwise polarized: 



/ (ExB) = +e3^^2^osV, 



(ExB) = +e3 M^A^cosV 



(10) 



Two waves are linearly independent and orthogonal to each other: 

, B^B^^ = . 



In the same manner we can solve a more general problem of constructing plane wave solutions 
with arbitrary wave vector k. Let us start with a scalar wave 



^ = A sin(A;oXo — kiXi) = A cos ip . 
The matrix of solutions Sbh will take the form 



(11) 



iko 


-ki 


-k2 


-k3 


h 


iko 


+k3 


-k2 




-k3 


iko 


+ki 


k3 


+k2 


-ki 


iko 



cos (f . 



(12) 



We have four formal solutions \E'", but they cannot be regarded as physical because each of them 
has a non- vanishing zero-component. However, we can use linearity of the Maxwell equation and 
combine elementary columns in (jl2p with any coefficients. In this way we are able to construct 
physical solutions. For shortness let us omit the factor A simp and operate only with the 
columns of the matrix 



iko 


-ki 


-k2 


-k3 


ki 


iko 


+k3 


-k2 


k2 


-k3 


iko 


+ki 


k3 


+k2 


-ki 


iko 
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Taking into account properties of the wave vector 

fc2 _ k2 = =^ k = fcon , = 1 , 
previous matrix can be rewritten as follows (the common factor /cq is omitted) 



(13) 



I -ni -n2 -ns 

ni i +713 —"2 

n2 — na i +ni 

na +712 —ni i 



I —ni 

711 i 

712 — J^a 



-7i2 -ns 

+7l3 -712 

^ +ni 



na +7i2 — ni 



First, with the help of the column (0) let us produce a zero at the first component of the columns 
(1) - (2) - (3): 



i 

ni —iniui + i 

712 — mi7l2 - ^3 

ns — mi7i3 + 712 











-in2ni + 77,3 
-m27i2 + i 
-m27i3 — ni 



-m^ni — 77,2 
-m37i2 + ni 
-in^n^ + i 



Noting that the column (3) is a linear combination of the columns (1) and (2): 

i 712 (1) - i ni (2) = (3) ; 

in other words, solution (3) is a linear combination of (1) and (2). Therefore from the above it 
follows (multiplying the columns (1) and (2) by imaginary i): 



i 

7ii n\ — 1 n2ni + i7i3 

7Z2 7Xi7l2 — m3 7X27X2 — 1 

713 nin^ + m2 7i2«3 — ini 



(14) 



Thus, two physical Maxwell solutions are 







1 

77,177,2 - in^ 

riius + in2 





+ icB^ 
E^ + ciB^ 





7727X1 + 7713 

7127X2 - 1 
7X27X3 - XTXl 





E^ + icB^ 
E^ + icB'^ 
E^ + ciB^ 



(15) 



or differently (the factor A cos (p is omitted) 

I E = (7xf — 1, 7X17X2, 7X17X3) . 

II E = (7X27X1, 7X27X2 — 1, 7X27X3) 



cB = (0, -7x3, 7x2) ; 
cB = (7x3, 0, -7x1) . 



It is the matter of simple calculation to verify the identity for amplitudes: cB 
these waves are clockwise polarized: 



(16) 

E. Also, both 









ei 62 


63 








E^ X B^ = 




n\ 


- 1 7117X2 


7X17X3 


= (1 


- n\) (7x161 + 7x262 + 7x361) 


~ k , 








-7x3 


n2 








E^^ X B^^ = 




ei 


62 


63 








7X27X1 


7X27X2 — 1 


7x27x3 


= (1 


- 7x2) (7x161 + 7x262 + 7x361) 


~ k 






7x3 





-7X1 









Besides, they are independent solutions (but not orthogonal ones) 

E-^E-^-^ = -7X17X2 , B-^B^^ = -7X17X2 • 



(17) 



(18) 
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5 Dual symmetry of Maxwell equations 

Let us consider the known dual symmetry in matrix formalism: 





E + icB 



. 



It is evident that there exists a simple transform, multiplication by imaginary i, with the fol- 
lowing properties: 



-ido + a^dj) = , 











+ icB^ 




iE - cB 



E 



D 



-cB , cB^ = +E . 



(19) 



It is the dual transformation of the electromagnetic field. Several points should be clarified. 
First, this transformation is not a symmetry operation in presence of external sources. In this 
case we have 



{-ido + a^dj] 



and further 



= i * , 


1- 


= -i^^ , 










E^ + icB^ 




iB - cB 






1 


P 






B + icB 


eo 








{-ido 


+ dj 




1 


ip 






eo 


-j 


E^ 


= -cB , 


cB^ 


= +E 



(20) 



(21) 



Second, to save the situation one can extend the Maxwell equations by introducing magnetic 
sources: 






1 


Pe + ipm 


E + icB 




y e ~l~ Jm 



{-ido + oPdj] 

which permits us to consider the dual transformation as a symmetry: 



(22) 



= i^ , ^ = _i , 



-ido + a^dj) = — 



1 

eo 










+ icB^ 




K - cB 






■ n 




Prm Je 



E 



D 



Pm ~l~ ^Pe 
y m Je 

-cB , cB^ = +E 



D 



+Jm 5 Pm — ~^Pi 



3m Je 



(23) 



In real form eqs. (I23p will look 



do — — \- oi-'dj- 



d 



2i 
^ + ^ 



'0- 

Re{J) = — 
eo 



2 +«'^^- 2^ 



^ 2 2 
^' - J - J* 



Pe 
Jm 



Im{J) = — 
eo 



2i 

Pm 
je 
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that is 

doB + a^djE = Re (J) , d^E + a^djB = Im (J) . (24) 
Eqs. (|24p in vector notation read 

div E = — , rot E = H — , 

Co eo Oct 

div cB = ^ , rot cB = + 1^ . (25) 
eo eo Oct 

Let us turn again to the Maxweh equation without sources and examine the action of this 
transformation on a plane electromagnetic wave along the axis z (see ©): 

/co 

I E"^ = —k^A cos(A;o2;o — ^32^3) , = — ^ cos(A;oXo — ^3X3) . 

c 

After the dual transformation it becomes 

Eo = -koA cos(A;o2;o - ^3X3) , cBI, = -k^A cos(A;oXo - A;32;3) , 
It should be noted that the dual solution coincides with the wave of the type II according to ([9]) : 

II E^ = k^A cos(A:oXo — k^x^) , B^ = —A cos{koxo — k^x^) . 

c 

In other words, the dual symmetry provides us with the possibility to construct a new linearly 
independent solution on the base of the known one. 

One addition should be made: the well-known continuous dual symmetry looks as a phase 
transformation over complex variables: 

e'^(E + icB), e-'^ (ije+jm) , e-'^ {pe + i Pm) ■ (26) 



6 On separating physical solutions of the Maxwell equations 
(real-valued scalar function $ ) 



Let us consider four formal solutions of the Maxwell equations 



iFo 


El 


F2 


F3 


-Fi 


iFo 


-F3 


F2 


-F2 


F3 


iFo 




-F3 


-F2 


Fi 


iFo 



Fa{x) = da^{x) 



(27) 



Physical solutions should be associated with the following structure: 



E + icB 

when zero-component of the 4x1 column vanishes. 

Let the function <I>(x) be taken as real- valued. We should find all possible solutions to the 
following equation: 





E + icB 



(28) 
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Let us separate real and imaginary parts in Ac : Ac = flc + ^^c- Relation (j28|) takes the form 

(ao + ibo)iFo + (ai + ibi)Fi + (03 + ^62)^^2 + (^3 + i&3)i^3 = , 
-(ao + ibo)Fi + (ai + ibi)iFo - (02 + i62)i^3 + {as + ibz)F2 = Ei + icBi , 
-(ao + ibo)F2 + (oi + ibi)F3 + (02 + i&2)«i^o - (^3 + ^^3)^1 = E2 + icB2 , 
-{ao + ibo)F3 - {ai + ibi)F2 + (02 + ib2)Fi + (03 + ib3)iFo = E3 + icB^ , 

from whence it follows 

- 60F0 + aiFi + 02^2 + asFg = , agFo + biFi + 62F2 + 63F3 = , 

-ao-Fi - biFo - 02^3 + 03-^2 = El , -boFi + aiFo - 62-P3 + ^3-?^2 = cBi , 

-ao-F2 + aiF^ - 62-^0 - 03^1 = E2 , -boF2 + 61 F3 + 02^0 - 63F1 = ci?2 , 

-00^3 - aiF2 + 02^1 - 63F0 = , -60F3 - 61F2 + 62F1 + 03^0 = cBs . (29) 

Let us consider equation rot E = —OqcB. Taking two identities 

51^2 - d2Ei = di[-aoF2 + aiF^ - 62-^0 - a-sFi] - 
-d2[-aoFi - biFo - 02-^3 + 03-^2] = 
= aidiF^ - b2diFo - a^diFi + bid2Fo + 0202^3 - 0392-^2 , 
-O0CS3 = 60^0^3 + bidoF2 - b2doFi - a^doFo ; 

we produce an equation 

aidiF^ - 62^1^0 - a^diFi + 6i92-Fo + a2d2Fs - 03^2^2 = 
= bodoFs + bidoF2 - b2doFi - asdoFo , 

from whence substituting identity OoFq = diFi + c?2-F2 + d^F^, we obtain 

aidiFs + 0292^3 = bodoFs - 03^3^3 ; 

that is 

[ bodo - {aidi + 02^2 + 03^3) ] F3 = . (30) 
In the same manner we get 

[ bodo - {aidi + 02^2 + 03^3) ] Fi = , 

[ bodo - {aidi + a2d2 + 0353) ] F2 = . (31) 
Now consider equation rot cB = 9oE. Calculating two rerms 

dicB2 - d2cBi = 

= bidiF^ + a2diFo - b^diFi - aid2Fo + 62^2^3 - 63^2^2 , 
doE-i = -aodoFs - 01^0^2 + a2doFi - bsdoFo = 
= -aodoFs - aidoF2 + a2<9oFi - 63^1^1 - 63'92F2 - 63^3^3 , 



we arrive at 



bidiFs + a2diFo - b^diFi - aid2Fo + 62^2^3 - 63^2^2 = 
-aodoFs - aidoF2 + a2«9oFi - 6391^1 - 63^2^2 - 6393F3 , 
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or 

[ aodo + ihdi + 62^2 + 63^3) ] F3 = . (32) 

Analogously, we get 

[ aodo + ihdi + 62^2 + hds) ] Fi = , 

[ aodo + ihdi + 62^2 + 63^3) ] F2 = . (33) 

Now let us consider equation div E = 0: 

= diEi + d2E2 + dsEs = -aodiFi - hdiFo - a2(9iF3 + a3(9iF2 - 
-00^2^2 + aid2F3 - b2d2Fo - 03^2^1 - 0093^3 - 0153^2 + 0293-^1 - 63^3^0 , 

that is 

- aoidiFi + ^2^2 + c>3F3) - ihdi + 62^2 + 63^3)^0 = , (34) 
which is equivalent to 

[ aodo + ihdi + b2d2 + 63^3) ]Fo = 0. (35) 
It remains to consider equation div B = 0: 

= -bodiFi + aidiFo - hdiF'i + 6391F2 - 6o92-F2 + ^i92-F3 + a292-Fo - ^3^2 -Fi - 

-bodsFs - 6193F2 + b2dsFi + a^d^Fo , 

or 

- boidiFi + ^2^2 + d^Fs) + iaidi + 02^2 + 0393)^0 = , 
which is equivalent to 

[60^0 - iaidi + 02^2 + 03^3)] Fo = . (36) 

Thus, to construct physical solutions of the Maxwell equations as linear combinations from 
non-physical ones (see (|28|) ) 



(ao + i6o)^° + (ai + ^^i)^^ + (a2 + ifc2)*^ + (as + ^^s)^^ 





E + zcB 



(37) 



one must satisfy the following 8 equations: 

[ bodo - iaidi + 02^2 + 03^3) ] Fc = , 

[ aodo + ibidi + 6252 + bsds) ] F, = ; (38) 
where c = 0, 1, 2, 3; F^ = dc<^. 
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7 On separating physical solutions of the Maxwell equations 
(complex-valued scalar function $ ) 

Let us consider four formal solutions of the Maxwell equations 





Fi 


F2 


F3 


-Fi 


iFo 




F2 


-F2 


F3 


iFo 


-Fi 


-F3 


-F2 


Fi 


iFo 



Faix)=da^x) . 



(39) 



Let the function <I>(x) be taken as complex- valued. We should examine relationship for Xa 
defining all possible solutions to the following equation: 





B + icB 



Let us separate real and imaginary parts in Ac and <I>(x) and Fc{x): 

Xc = ac + ibc, ^(x) = L{x) + i K{x) , Fc{x) = Lc{x) + i Kc{x) 
Relation ()40p takes the form 

(ao + ibo)i{Lo + iKo) + (ai + ibi){Li + iKi) + 
+(02 + ib2){L2 + iK2) + (as + ib3){L3 + iKs) = , 

-(ao + ibo){Li + iKi) + (oi + ibi)i{Lo + iKo) - 
-{(12 + ib2){L3 + iKs) + {as + ib3){L2 + iKi) = Ei + icBi , 

-(ao + ibo){L2 + iKi) + (ai + ibi){L3 + iK^) + 
+{a2 + ib2)i{Lo + iKo) - (ag + ib3){Li + iKi) = E2 + icB2 , 

-(ao + ibo){L3 + iKg) - (ai + ibi){L2 + ii^a) + 
+(02 + ib2){Li + iKi) + (ag + ib3)i{Lo + iKo) = E3 + icBs , 

from whence it follows 



: 
: 

El - 
Bi ■■ 

E2- 
B2 

E3- 
B3 



-boLo + aiLi + a2L2 + asLa) + (-ao-f^o - biKi - b2K2 - bsKs) 
+aoLo + biLi + 62-^^2 + bsLs) + (-60-K'o + aiKi + a2i^2 + o-a-^^s) 

-aoLi - biLo - a2L3 + a3L2) + {+boKi - aiKo + 62/^3 - ^3-^2) 
-60^1 + aiLo - b2L3 + 63L2) + {-aoKi - biKo - a2i^3 + a3i^2) 

-aoL2 - 62-^-0 + aiL3 - a^Li) + (+60-^2 - a2-f^o - biK^ + b^Ki) 
-boL2 + a2Lo + biLs - 63L1) + (-aoi^2 - &2i^o + - 03-^^1) 

-aoLs - 63L0 - aiL2 + a2Li) + (+6o-?^3 - o-a-f^o + biK2 - b2Ki) 
-60L3 + a3Lo - biL2 + b2Li) + (-ao-f^s - bsKo - aiK2 + a2i^i) 



(40) 



(41) 



Substituting these expression into Maxwell equations and performing calculation like in 
previous section we get: 
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div E = 


=^ 


-{aodo + bjdj) Lq + (feo^o - ajdj) Kq 


= 0, 




div B = U 




+(00(70 - ajOj) Lq + (aodo + bjdj) Kq 


= , 


92^3 


- ^3^2 = +doEi 


=^ 


-{aodo + bjdj) Li + {bodo ~ ajdj) Ki 


= 0, 


92-^3 


— dzE2 = —BqEi 




+{bodo - ajdj) Li + (ao^o + bjdj) Ki 


= 0, 




- diBs = +doE2 




-{aodo + bjdj) L2 + (60^ - ajdj) K2 


= 0, 




— diE^ = —OqB^ 




+(60^0 - ajdj) L2 + (ao^o + bjdj) K2 


= 0, 


diB2 


- d^Bi = +doEs 




-{aodo + bjdj) L3 + (6o9o - ajdj) 


= 0, 


diE2 


— 52-Ei = —doBs 




+{bodo - ajdj) L3 + (ao^o + bjdj) Ks 


= 0, 



Thus, for physical Maxwell solutions the following equations must hold: 
- (aodo + bjdj) Lc + {bodo - ajdj) Kc = ; 

+{bodo - ajdj) Le + (ao^o + bjdj) Kc = 0. (43) 

In particular, we have two more simple equations when taking real or imaginary scalar 
functions: 

^ = L + iO , 

ac, be <^ (ao i9o + bj dj) Lc = , {bo do - aj dj) Lc = 
^ = + iK 

a'cX ^ {a'odo + b'jdj)Kc = 0, {b'o do - a'j dj) = . 

(44) 

In the simplest case of a plane scalar wave 

$ = e*('=o^°-*^'^') = cos + i sin 99 , 

ip = ko - P = kcx" , (45) 
Lc + iKc = —kc sin if + i kc cos 

previous equations take the form 



(ao ko + bj kj) kc = , {bo ko — aj kj) kc = 
{a'o ko + bj kj) fee = , {b'o ko — a'j kj) fee = . 

(46) 
that is 

a'c = ac , He = be , 
{ao ko + bj kj) kc = , {bo ko — aj kj) kc = (47) 



$ = L + z , 

ac, be 
^ = + iK 



a'cX 
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In general case equations for Oc, be and a'^, b'^ may not coincide. 
Let us turn again to eqs. (|44p and translate them to variables 

jp* -\- Fc jp* — 

c 2 2 

^* _l_ ^ p* — p 

- (ao^o + ftjSj) -^^ — - + {bodo - ajdj) i — - = ; 

+{bodo - ajdj) El±I± + (ao^o + 5.9.) i ^^^^ = . 

or 

1 1 

- [ -{aodo + bjdj) - i {bodo - ajdj) ] Fc + -[ -{aodo + bjdj) + i {bodo - ajdj) ] F* 
^ [ +{bodo - ajdj) - i {aodo + bjdj) ] Fc + ^[ +{bodo - ajdj) + i {aodo + bjdj) ] F* -- 
or 

\ [ -{aodo + b,dj) - ^ {bodo - a,d,) ] + i [ -{aodo + 6,9,) + . (60^0 - a,d,) ] F^ 

- [ +i {bodo - ajdj) + (oo^o + bjdj) ] Fc + {bodo - ajdj) - {aodo + bjdj) ] F* -- 

Summing and subtracting two relations, we arrive at 

[ -{aodo + bjdj) + i {bodo - ajdj) ] F* = 
[ -{aodo + bjdj) - i {bodo - ajdj) ] = 

They can be rewritten as follows: 

[ - (ao + ibo) do + i{aj + ibj) dj]Fc = Q , 
[ - (ao - ibo) do - i{aj - ibj) dj ] P* = , 

or shorter 

i-Xodo + i Xj d,]F, = 0, 
[-Xod*o-t A* dj ]F* = 0, 

8 Separating physical solutions of the plane wave type 

Let us apply the general relations ([38]) to the case when 

^ = A sin 93 , if = koxo — /C3X3 , 
Fo = koA cos ip , Fi = , F2 = , F3 = -k^A cos ip . 

Eqs. (pSl) then give 

60^0 + 03^:3 = , aoko - 63/C3 = . 
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For a wave spreading in the positive direction = +/co > 0, and eqs. (j^U]) give 

bo = -as , 63 = ao , 
and correspondingly relationship (j37p looks 



E + icB 



(ao - iag)^'" + (03 + iao)^^ + (^1 + i^i)^'^ + (02 + ^^2)*^ 



(51) 



coefficients at ^'^ and ^'^ are arbitrary. To understand this fact let us recall the explicit form of 
- see (El): 



cos (f . 



iko 








-ko 





iko 


ko 








-ko 


iko 





ko 











One may separate two subsets of non-physical solutions : 

(ao - iasj^-" + (03 + iao)^'^ = (ao - ia3)A:o < 
so relationship ([5T]) reduces to 

(ai + ibi)^^ + (a2 + ^62)^^ = (ai + i6i 



i 




-1 













+ i 





+1 




i 



(52) 





iko 

-h 





(a2 + ib2) 





iA;3 




(53) 



Now let us consider a more general case when 

$(x) = Asin(A;oXo — kx) = A sin 99 , 
Fo = koAcosip , Fi = —kiAcosip , 
F2 = —k2A cos if , F3 = —ksA cos 99 . 

Eqs. §Si) give 

60^0 + ai^i + 02/1:2 + aa/ca = , aoko - biki - b2k2 - b^k^^ = 



(54) 
(55) 



and additionally the identity ko = + ^2 + ^3 holds. One may introduce parametrization 

kj = koUj, rijUj = 1, then eqs. ([55]) read 



bo = -(aini + 02722 + a^ns) , ao = {biui + 62^2 + fts^-s) • 

Now turning to ([37]) and excluding variables ao, 60 one gets 

[(6ini + 627^2 + bsn^) - i{aini + a2n2 + 03713)]^'° + 
+(ai + ibi)^^ + (a2 + ^^2)^^ + (03 + ^^s)^^ = ^ • 

Taking into account (see ([2] 



/C0^COS(/3 , 



(56) 



(57) 





-ni 


-712 


-773 


77-1 


i 




-772 


77-2 




i 


771 






-ni 


i 
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from ([57|) we get (factor kQA cos ip is omitted) 



[{bi - iai)ni + (62 - ia2)n2 + (63 - ^03)713] i - {ai + ibi)ni - (02 + 162)^^2 - (as + ^&3)"'3 
[(61 - iai)ni + (62 - ia2)n2 + (63 - ^a3)"-3] "-i + (oi + + (^2 + i&2)"-3 - (^3 + i&3)"'2 
[(61 - iai)ni + (62 - ia2)n2 + (63 - «a3)"-3] "-2 - (ai + i6i)n3 + (02 + ib2)i + (03 + «&3)"'i 
[(61 - iai)ni + (62 - ia2)n2 + (63 - ^03)123] 713 + (oi + ibi)n2 - (02 + 762)^1 + (03 + 763)7 



and further 



[(61 - 701)771 + (62 - 702)772 + (63 - ^03)773] 
[(61 - 701)771 + (62 - ia2)n2 + (63 - ia3)n3] 
[(61 - 701)771 + (62 - ia2)n2 + (63 - ia3)n3] 





ni + (oi + 761)7 + (02 + 762)773 - (03 + 763)772 

772 - (oi + 761)773 + (02 + 762)7 + (03 + 763)771 

ns + (ai + 761)772 - (02 + 762)771 + (03 + 763)7 



Let us introduce three elementary solutions, associated with coefficients (uj + 76j): 



(1) 



(2) 



(3) 





{bl- 
(61- 
(61- 



(62- 
(62- 
(62- 



(63 

(^3 

(63 



70l) 77l77i + (oi + 761)7 
70l) 771772 - (Oi + 761)773 
70l) 771773 + (Oi + 761)772 



7O2) 772771 + (02 + 762)773 
7O2) 772772 + (02 + 762)7 
7O2) 772773 - (02 + 762)771 



703) 773771 - (03 + 763)772 
703) 773772 + (03+763)771 

^03) ns 773 + (03 + 763)7 



,E(i) 



,E(2) 



,E(3) 



61771771 — 61 

61 771772 — 01773 

61 771773 + Ol772 



62772771 + O2773 

62772772 - 62 
62^2773 - 02711 

6377.3771 - O3772 
&3"-3f^2 + 03"-! 
&3"-3 "-3 - ^3 



,cB(i) 



, CB(2) 



, CB(3) 



-Oi77i77i + Oi 
-Oi77i772 - 6in3 
-Oi77i773 + 6in2 

-02772771 + 62773 
-02772772 + O2 
-02772773 - 62711 



-O3773771 - 63772 
-O3773772 + 63771 
-O3773 773 + 03 



Let us show that three types of solutions are linearly dependent. It suffices to examine their 
linear combinations (for definite consider electric field): 



Ai E(i) + A2 E(2) + A3 E(3) = , 



that is 





6i77i77i — 61 




62772711 + O2773 




63773771 - 


- 03772 









61 771772 - O1773 


+ A2 


62772772 - 62 


+ A3 


63773772 - 


f O3771 









61 771773 + 01^2 




^2?^2"-3 - a2ni 




63773 773 


-63 








It remains to show that the determinant of the 3x3 matrix vanishes 
6i77i77i - 61 b2n2ni + 02773 63773771 - 03712 

det 6i77i772 - 01^3 62772772 - 62 63/73772 + 03^1 = 

6i77l773 + 01^2 62772773 - 02771 63773 773 " h 



(60) 
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Let biui = si, 62"'2 = ■52, b^n^ = S3 then eq. ([5n|) looks 



sini — bi 
sin2 — ains 
Sins + ain2 



S2ni + 02723 
5271-2 - 62 

S27J3 — a2ni 



53711 - 03712 

53712 + 03711 

S3"'3 - ^3 



= (5irii - 6i)[(s27i2 - &2)(s3n3 - fts) - (s3"'2 + a^rii){s2n'i - 02711)] - 

-(5i7l2 - Oin3)[(52ni + 02n3)(s37l3 - 63) - (ssTli - 037l2)(52ri3 - 027li)] + 
+ (5i7Z3 + Oi7l2)[(5271l + 027l3)(53n2 + 03711) - (53711 - 037l2)(52n2 " ^2)] = 



= (5ini - 6l)[5253?T.2n3 - 527X263 - 537X362 + 6263 - S2S3".27).3 + 537Xi7l202 - 7X17130352 + 02037X?] - 
-(5i7X2 - 011X3) [52537X11X3 - 527X163 + 02537X3 - 027X363 - 52537X17X3 + 53027X1 + 03527X27X3 - 03027X17X2] + 
+ (S17X3 + 017X2) [52537X17X2 + 52037xf + O2537X27X3 + O2O37X17X3 - 52537Xi7X2 + 53627X1 + 03527X3 - O3627X2] 



= -517X1527X263 - 517X1537X362 + 517X16263 + 5l7Xi537Xi7X202 - 5l7Xi7Xi7X30352 + 5l7Xi02037xf - 
+61527X263 + 61537X362 - 616263 - 61 537X11X202 + 617X17X30352 - 6i02037X^ + 

+ 5l7X2527Xl63 - 5in20253n| + 517X2O27X363 + 5l7X253027Xi - 5 1 7X2 O3 5 2 7X2 7X3 + 5l7X203027Xl7X2 - 
-O1IX352IX163 + Oi7X302537x| - O17X3O2IX363 + O17X353O27X1 + O1IX3O3527X27X3 - O17X3O3O27X17X2 + 

+517X352037X1 + 51X1302537X27X3 + 517X302037X17X3 + 511X353627X1 + 51X1303527X2 - 517X303627X2 + 

+oi7X252037Xi + 017x202537x2x13 + 017x20203x117x3 + 017X253627X1 + 01X1203527X2 - ain2a^b2n2 
and further 

= -6162637x17x3 + 616263x^1 + 616302x117x2x13 - 6162037x1x127x3 + 61O2O3XX1 - 
+6162637X2 + 616263XXI - 616263 - 616302XX17X2XX3 + 6162037X17X2XX3 - 610203XX1 + 

-6163O2XX1XX2IX3 + 6i6302nin2XX3 + 6i63a2XX2XX3XXi - 6162O3XX1XX2XX2XX3 + 61O3O2XX1IX2 - 

-626301X137X2X11 + O163O2XX3 - 016302X13 + 63O1O2XX37X1 + 62O1O37X2XX3 - Oi02a3XXiXl2XX3 + 

+6162037x21x37x1 + 6ia263XX3XXixx2 + 61O2O37X1XX3 + 616263XX3XX1 + 6162O3XX1XX3XX2 - 6103627x1x127x3 + 

/I O 

+0162O37X2XX1 + 01O263XX2XX3 + 0102037x17x27x3 + 0162637x17x27x3 + 01O3627X2 - 01O362IX2 == ; 

it is easily verified that all terms cancel out each other indeed. 

Let let us consider one other example and start with a complex scalar plane wave: 

^x) = e^^'^o^"-'^^^:') , Lc = kc cos ^ , Kc = kc sin 99 ; 

and eqs. (I43p take the form 

— (oo^o + bjkj) kc cos if + (6oA;o — o.jkj) kc sin 99 = ; 
+(6oA;o — cLjkj) kccoscp + (oq/co + bjkj) kcsiiup = . 

from whence it follow 

boko + (oi/ci + 02^)2 + 03/^3) = , oo^o - (6iA:i + 62/02 + 63/^3 = . 

or 

60 = -(oixxi + 02XX2 + 03XX3) , 00 = +(6ixii + 62x12 + 63XX3) . (61) 
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Correspondingly, eq. ([57|) gives 

^{x) = [+{bini + 62712 + hn^i) - i{aini + 02^2 + 03713)] + 
+(ai + ibi)^'^ + (02 + ib2)^'^ + (03 + ^63)^'^ 





E + icB 



where (see (E 



Further we get 



« —Til 


-712 


-713 


Til i 


713 


-772 


712 -713 




771 


^3 712 


-Til 


i 



ikoA (cos + 7 sin 99), 



^^) = [-i(a + 7b) n + (ai + 761)^'^ + (02 + 762)^^ + (^3 + «&3)^^] 
(a + 7b)n - (oi + ibi)ni - (02 + 762)?^2 - (03 + 763)713 
-7(a + 7b)n 771 + (ai + ibi)i + (02 + 762)713 - (03 + 163)712 
-i(a + 7b)n 772 - (ai + 761)713 + (02 + 162)1 + (03 + 163)711 
-i(a + ib)n 113 + (ai + i6i)?i2 - (02 + 162)111 + (03 + ^^3)^ 



ikoA (cos ip + i sin ip), 



that is 





(a + ib)nni - (oi + 161) + i(a2 + 162)113 - 1(03 + 163)112 
(a + ib)nn2 - i(ai + 161)113 - (02 + 162) + 1(03 + 163)111 
(a + ib)nn3 + i(ai + 161)112 - i(a2 + 162)111 - (03 + 163) 



koA (cos ip + i sin ip) 



which is equivalent to 







(aniii - ai - 62113 + 63112) + i(bnni - 61 - 11203 + "-3^2 
(ann2 - 02 - 63111 + 61113) + i(bnn2 - 62 - 11301 + nia2 
(ann3 - 03 - 61112 + 62111) + i(bnii3 - 63 - nia2 + n2ai 



k^A (cos 93 + i sin ip) 



C = n (nb) — b + a X n, 



With the use of notation 

L = n (na) — a — b x n , 
relationship (j62p can be written shorter 

E + icB = (L + iC) koA (cos (p + i sin (p) , 

from whence it follow 

E = koA (cos LpJj — sin cpC) , cB = k^A (sin (/?L + cos (pC) 

One can readily prove identities: 

l2 = = a^ + b^ - (na)2 - (nb)^ + 2n (a x b) , 
LC = 0, EB = 0, e2 = c^B^ , 
Ln = 0, Cn = 0,En = 0, Bn = . 
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General expressions for L, C may be decomposed into the sum: 

L = Li + L2 , C Ci + C2 , 
/ b = : Li = n (na) — a , Ci = a x n , 

El X cBi = k^A^ (Li X Ci) = klA^ [ - (na)^ ] n ; 



// a = : L2 = -b X n , C2 = n (nb) - b , 

E2 X cB2 = klA^ (L2 X C2) = klA^ [ - (nb)2 ] n , 

In other words, these two electromagnetic wave are clockwise polarized. 
For a particular case when b = a, we get 

I Li = n (na) — a , Ci = a x n , 
El = koA{cos(p Li — simp Ci) , 
cBi = koA (sint^Li + cos(^Ci) ; 
II L2 = —a X n = — Ci , C2 = n (na) — a = Li , 

E2 = koA{— sin iphi — cos (p Ci) , 
CB2 = koA (— sin tpCi + cos tplii) . 

so constructed waves are linearly independent and orthogonal: 

E1E2 = , B1B2 = . 



(66) 



(67) 



(68) 



In view of linearity of the Maxwell equations any linear combination of the type is a solution as 
well: 



E — ci El + C2 E2 



B = ci Bi + C2 B2 



(69) 



9 Cylindrical waves 

Let us start with a cylindrical scalar wave (below for brevity ko = E): 



^lExo „ikz im4 



R{p), 



(70) 



Xi = p cos (f) , X2 = p sin (f) , xs = z , 
d sin ^ 5 



d _ dp d ^ dcj) d _ 

dxi dxi dp 8x2 d(p dp p dcp ^ 

d dp d ^ dcp d ■ ^ ^ ^ 

dx2 dx2 dp dx2 d(f) dp p dcf) 



(71) 



Corresponding electromagnetic solutions are to be constructed on the base of relations: 





E + icB 



(72) 



Let us find Fa- 



iE e^-^^o e^*^^ e'™"^ R{p) 



F3 = ik e^^^o e'^"" e'""^ R{p) , 



(cos (p im ) H , 

dp p 

, . , d . cosd). „ 

(sm m — +im ) it , 

dp p 



(73) 
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The main equations to solve are 

(-Ao^o + iXjdj) Fe = (74) 

or 

[-i\oE-\^k + i\x{cosct)- ^ +^A2 sin0 — + T i^;^ = q , 

op P ocp op p 0(p 



or 



r , „ , , , ^ . ,s d — Ai sin q!) + A2 cos o!) , „ 

-zAoF-A3A; + z(Aicos.^ + A2sin(^ — +i ^ r ]Fc = 0. 75 

op p o<p 



Let c = 1: 



Let c = 2: 



5 

[— iAo-E — A3A; + i(Xi cos ^ + A2 sin (f)) - — h 

op 

p ocp dp p 



d 

[—iXqE — X^k + iiXi cos + A2 sin (f)) — — h 

Op 

-Aisin^ + A2cos0 9 j^j, d cos(f) 

+t zi e«™-<P (sin0— + tm -)R = 0, 

p Ocp dp p 



Noting that 



dp p 

Fi - zF2 = e^-^^° e''^^ e'^rn-i),p^d ^ m ^ 

dp p 

from previous two equations we get 

, „ > , , , ,s d Ai sin cb — Xo cos , ^ , d m. ^ 

[-iXqE - Xsk + i Ai cos </> + A2 sin (/-^ + m + 1 )R = 0, 

op p dp p 

r X „ X , , X . , N 9 Ai sin ch — X2 cos ch , ^ , , d m , „ , , 

-^AqF - A3A; + i{Xi cos <^ + A2 sin <^ — + ^ {m-l)]i— + —)R = 0. 77 

op p dp p 

In turn, when c = 0, 3, we get one the same equation: 

. , „ , , , . . ,^ 9 Ai sin (i) — A2 cos c!) , , 

[ -iXoE - A3A; + i{Xi cos ^ + A2 sin ^ + ^ ^ m ] i?(p) = , (78) 

We readily note that there exist very simple way to satisfy these three equations on param- 
eters Ac: 

- iXoE - A3A; = , Ai = , A2 = . (79) 
Let us demonstrate that no other solutions exist. To this end, with notation 

C = —iXqE — X^k , ^ = Ai cos i;^ + A2 sin cj) , B = Xi sin cp — X2 cos cp . 
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let us rewrite eqs. ([77|) - ([75|) in the form 



^ ■ A d B , . ^ .dR m 

C + lA— + - m + 1) (- R) = Q 

dp p dp p 

^ ■ A d B . ^ .dR rn ^. 

C + lA— + - m - 1 (— + — i?) = 
dp p dp p 



d B 

[C + iA— + —m]R = 0. (80) 
dp p 



Combining two first equations in ([80]) (summing and subtracting), we get 

dR d'^R niBdR mB 

C— + iA -— + ^R = , 

dp dp-^ p dp p^ 

mC ^ , d R BdR ^B ^ 

R - imA \ ; m^—R = , 

p dp p p dp p^ 



^ „ . dR B ^ , . 

CR + iA — + m — R = 0. 81) 
dp p 

After differentiating the third equation will look 

_ dR , d'^R mB dR mB 

C — + iA—y + 1 ^ = 

dp dp'^ p dp 

which coincides with the first equation in (ISip . Therefore, the system ()8ip is equivalent to 

mC „ d R BdR oB ^ 

R - imA \ ; m^^R = , 

p dp p p dp 

C R + iA*^ + m — R = Q . (82) 
dp p 



The system (??) can be rewritten as follows: 



C B dR mA m , dR B 

— R+ — + i^rR {iA — + m — 

p p dp p'^ p dp p 



CR + iA^ + m — R = 
dp p 



here the first equation with the use of the second equation gives 

B dR mA 

- — + i^R = 0. 

p dp p'^ 



Therefore system (|82p is equivalent to 



dR mA 

B — + i = , 

dp P 

iA'^ + m- R + C R = 
dp P 



which in turn is equivalent to 



d 777 

{B + iA)i^ + -)R + C R = 0, 
dp p 

iB-iA){^-—)R-CR = 0, (83) 
dp p 
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Noting identity {A + iB){A - iB) = Af + A| , one reduces eqs. (j83|) to the form 



(A? + A^)(— + -)R +{B- iA)C R = 0, 
dp p 

n IT) 

{\l + XDij^ - -)R + {B + iA)C R = 0. 

Remembering that 

C = —iXoE — A3A; , A = \i cos (f) + X2 sin (p , B = Xi sin cj) — X2 cos (j) . 
we immediately conclude that eqs. ()84p can be satisfied only by the following way: 



(84) 



(a? + a; 



2) ; 



c = o. 



(85) 



In other words, This means that relations (|79p provide us with the only possible solution in 
terms of two (1x4) columns 



^> = Ao^'° + A3^'^ 



iFoix) 
-Fi{x) 
-F2{x) 



-XpE + iXsk = 
F3{x 

F2{X 

-Fi{x 
iFo{x 
sin 



d 



Fq = iE ^ , Fi = (cos <p im — 

dp p 



d cos (j) 

F2 = (sm (j) — — h im ) 4> 

dp p 



F^ = ik ^ 



Solution explicitly reads by components: 

(^')o = (-Ao ^ + i A3A;) $ = 



£3 + iBs = (^')3 = (-Ao ik - X3E) $ 
El + iBi = (^)i = 



-Ac 



E^-k^ 
E 



, , . d . smc/). ^ / . , • cos(/), 

Ao (cos im j + A3 (sm (p - — h im j 

dp p dp p 



E2 + iB2 = (^')2 = 
After simple rewriting we get 



. , . . , d . cos , , , , d . sm ( 

(— Ao (sm cp - — \- im j — A3 (cos cp im 

dp p dp p 



^3 + iBs = i^h 



-A3^^<J> 



El + iBi = ^ 

E2 +iB2 = ^ 



d sine 

-ik icosq) — rni 

dp p 



d cos ( 

+ E (sm (j) — — h im 

dp p 



, d . cos 4). ^ , , d . sm0, 

i—ik (sm (h - — h im j — E (cos (p im j 

dp p dp p 



For special cases we get most simplicity: 

k = +E , 



E3 + iB3 = 0, 



(86) 
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Ei+iBi = A3 
E2+iB2 = A3 



d smcp d coscj) 

—I (cos (j) im ) + ( sm (p - — h im ) 

dp p dp p 

d cos(f) d sin0 

[—1 (smcp - — h im ) — (cos im ) 

dp p dp p 



$ =-iA3 e^'^i^--)^, 
dp P 

$ = -A3 e^^(^-^)$ , ; 
dp p 



k = -E , 



Es + iB3 = 0, 



Ei + iBi = A3 
E2+iB2 = A3 



d sm (j) d cos (p 

+1 (cos 6 tm ) + (sm 6 - — h tm 

dp p dp p 



$ = iA3 e-'H^ + -)^ 
dp p 



d cos( 

i+t (sm0 - — h im 

dp P 



d sine/) 

(cos (p im j 

dp p 



d 



A3 e-'^{—+-) $ 
dp P 



It seems reasonable to expect further developments in this matrix based approach to Maxwell 
theory, as a possible base to explore general method to separate the variables for Maxwell 
equations in different coordinates. Also it would be desirable to extent this approach to Maxwell 
theory in curved space - time models. 



Authors are grateful to all participants of scientific seminar of Laboratory of Theoretical 
Physics of Institute of Physic on NASB for discussion and advice. 
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